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Abstract 

The energies and widths of several D° meson bound states for different nuclei are obtained 
using a D-meson selfenergy in the nuclear medium, which is evaluated in a selfconsistent 
manner using techniques of unitarized coupled-channel theory. The kernel of the meson- 
baryon interaction is based on a model that treats heavy pseudoscalar and heavy vector 
mesons on equal footing, as required by heavy quark symmetry. We find D° bound states 
in all studied nuclei, from 12 C up to 208 Pb. The inclusion of vector mesons is the keystone 
for obtaining an attractive D-nucleus interaction that leads to the existence of -D°-nucleus 
bound states, as compared to previous studies based on SU(4) flavor symmetry. In some 
cases, the half widths are smaller than the separation of the levels, what makes possible their 
experimental observation by means of a nuclear reaction. This can be of particular interest 
for the future PANDA@FAIR physics program. We also find a D + bound state in 12 C, but 
it is too broad and will have a significant overlap with the energies of the continuum. 



1 Introduction 



The D-meson nucleus optical potential has been a subject of intense study over the last years. In 
particular, modifications of the D meson properties in an hadronic environment might influence 
the rhythm of charmonium production [TJ [21 El SI El EJ [T] , as a complementary explanation for 
charmonium suppression in a Quark-Gluon Plasma [8] . Moreover, D-meson bound states in 208 Pb 
were predicted [H] relying upon a strong mass shift for D mesons in the nuclear medium based on 
a quark-meson coupling model [TO]. The experimental observation of these bound states might be 
however problematic since, even if they exist, their widths could be very large [11] as compared to 
the separation of the levels. In Ref . [11] , the D-meson potential was obtained using techniques of 
self-consistent unitarized coupled-channel theory adapted to the meson-baryon interaction [12j [13] , 
which followed and extended the works using chiral Lagrangians and the Lippmann-Schwinger 
equation initiated in Refs. [TU [15] . This work followed a scheme similar to those used in previous 
approaches on the spectral features of D mesons in symmetric nuclear matter pjjj [T7J [181 HH] • The 
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systematic inclusion of medium corrections to the scattering equation is crucial for the generation 
and modification of the A c (2595) resonance in the nuclear medium and, thus, eliminate the main 
source of uncertainty of earlier evaluations of the D nucleus optical potential. In fact, this has 
been also done for the 77-meson selfenergy in a nuclear medium in Ref. [2T| . using the vacuum 
amplitude of Ref. [22]. The 77 nucleus optical potential was evaluated in a selfconsistent way, as 
it was done for the antikaon case in Ref. [23]. Later this optical potential was used for studying 
the possible 77-nucleus bound states in Ref. [24] , where the energy dependence of the 77 selfenergy 
around the rjN threshold was taken into account and shown to be very much relevant for the 
problem. 

Here we will undertake a similar study for the D°— nucleus bound states using a recently 
generated D meson selfenergy in the nuclear medium |25j. This model incorporates heavy-quark 
symmetry in the charm sector improving in this respect with the recent t-channel vector meson- 
exchange approaches [HI [HJ [HI [20]. As a consequence, the pseudoscalar D meson and the D* 
meson, its vector partner, are treated on equal footing. This new scheme generates a broad 
spectrum of new resonant meson-baryon states in the charm one and strangeness zero [25] arid the 
exotic charm minus one [27J sectors. Furthermore, this framework allows to obtain simultaneously 
the properties of D and D* mesons in nuclear matter [25]. The in-medium calculation of Ref. [25] 
includes Pauli blocking effects on the nucleon and the D and D* selfenergies in a self-consistent 
manner. Moreover in this work, a novel renormalization scheme is introduced that guaranties that 
the nuclear medium corrections do not depend on the ultraviolet cutoff used to renormalize the 
free space amplitudes. Compared to previous results [TTj IT5] , the width of the D meson in nuclear 
matter is small with respect to the mass shift and, therefore, bound states for D mesons in nuclei 
might be expected. 

This paper is organized as follows. In Sect. [2] we introduce the P-nucleus potential from 
Ref. [22] • The results for the different P-nucleus bound states are discussed in Sect. [31 where we 
also compare with other microscopical models. The conclusions are drawn in Sect. [H 



2 The D-nucleus optical potential 

In Ref. [25] the selfenergy of the D meson is evaluated in nuclear matter at various densities, p, as 
a function of the D energy, q°, and its momentum, q, in the nuclear matter frame. It is calculated 
by means of 

iW,g;p) = / ^ [T^^(P^P) + 3T^^/ 2 )(P^P) 
jp<p F {iTiy 1 

where p and pp are the momentum of the nucleon and the Fermi momentum at nuclear density 
p, respectively. The quantity T^ N (P°, P) is the in-medium DN s-wave interaction, with total 
four-momentum (P°,P) in the nuclear matter frame, namely P° = q° + E^(p) and P = q+p. 
Here, isospin symmetry is assumed, and the amplitude is summed over nucleons up to the Fermi 
level. Since we are interested in finding bound states, we shall be concerned about the s-wave D 
selfenergy around the DN threshold. 

The in-medium interaction T^ N is obtained self-consistently including all coupled channels 
with the same quantum numbers. In the work of Ref. |25j, the Bethe-Salpeter equation is solved 
with sixteen coupled channels for 1 = 0, J = 1/2, twenty- two channels for 1 = 1, J = 1/2, 
eleven channels for 1 = 0, J = 3/2 and twenty for I = 1, J = 3/2. The medium effects 
in the scattering amplitude are the Pauli blocking on the intermediate nucleon states and the 
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Figure 1: Selfconsistent D selfenergy of Ref. [25] for zero momentum as a function of energy, and 
for four different densities. The solid lines stand for the fitted functions of Eqs. (j!])-©. The dashed 
line shows the selfenergy for p = p calculated in the low density limit (without selfconsistency) . 



selfenergy of the charmed mesons (D and D*) in the intermediate channels. The pion and baryon 
selfenergies are not considered. While the baryon dressing did not change the qualitative behaviour 
of the D meson in nuclear matter [IT], the coupling to intermediate states with pions is of minor 
importance for the DN and D*N dynamics in a dense environment [25]. Note the importance 
of the selfconsistency for Ur)(q , q;p), since the in-medium amplitude T^ N contains the D(D*)N 
channel, which depends in turn on H D r D *\(q°,q; p). 

The D selfenergy, scaled by 2m d, is displayed with points in Fig. [1] as a function of the D- 
meson energy around threshold. It is shown for various nuclear medium densities, p, and with the 
D meson momentum fixed to zero. 

The D selfenergy is evaluated in infinite nuclear matter. In finite nuclei we use the local 
density approximation (LDA), substituting p by p(r), which is the local density at each point in 
the nucleus taken from experiment. For the s-wave, as it is here the case, it was shown in Ref. [28] 
that LDA gave the same results as a direct finite nucleus calculation. Then, the LDA D selfenergy, 
Il£)(g ,r) = U D (q°, 0, p(r)) allows to define an energy-independent local optical potential, 

T/ / \ U D (q° = m D ,q = 0,p(r)) 

VD\r) — , [2) 

2m d 

from its threshold (q° = m^, q = 0) value. This prescription gives a potential of ( — 25 — il4) MeV 
at normal nuclear matter density po — 0-17 fm~ 3 as can be read off from Fig. [TJ This would mean 
that one can expect bound states with approximately —20 MeV binding and a half-width of about 
14 MeV. 

However, both the real and the imaginary parts of the D selfenergy, around the D-meson 
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mass, show a pronounced energy dependence, as can be appreciated in Fig. [TJ For instance, the 
real part at q° — mo = —20 MeV is about one forth of its value at q° = mo- Hence, a realistic 
determination of the D bound states should take this energy dependence into account. Thus, we 
use an energy dependent optical potential defined as: 

v (r rp\ np(q° = m D + E,q = 0, p(r)) 

where E = q° — is the D° energy excluding its mass. In order to use the potential defined 
above in Eq. Q, the results of Ref. [25] are parameterized in terms of analytical functions in the 
energy range —25 MeV < E < as (see the solid lines in Fig. [TJ, 

Re[U D (q° = m D + E,6 ; p)] = a(p) + b(p)E + c(p)E 2 + d(p)E 3 , 

lm[U D (q° = m D + E,6;p)} = e(p) + f (p) E + g (p) E 2 + h(p) E 3 , (4) 

with 

a(p) = (-84033.1 p/p - 25727.8 p 2 /Po + 14536.2 p 3 /Po) M eV 2 , 

b(p) = (-8654.6 p/po + 6475.6 p 2 /Po) MeV , 

c(p) = -341.2 p/p + 447.8 p 2 /p 2 , 

d(p) = (-5.648 p/po + 8.777 p 2 /p 2 ) MeV" 1 , 

e(p) = (-10616.5 p/po - 39532.6 p 2 /p 2 - 2489.6p 3 /Po) MeV 2 , 

f(p) = (-3356.6 p/po - 1337.44 p 2 /p 2 ) MeV, 

g{p) = -291.09 p/p + 136.656 p 2 /p 2 , 

h(p) = (-6.70 p/po + 5.20 p 2 /p 2 ) MeV" 1 . (5) 

In the next section we will solve the Schrodinger equation, with both the energy dependent [Eq. ([3])] 
and independent [Eq. (|2J)] potentials, to find bound states for different nuclei through the periodic 
table. Since the D-meson optical potential is much smaller than its mass, we expect the relativistic 
corrections to be tiny and certainly smaller than the theoretical uncertainties of the interaction. 
We will also discuss the implications of our results in the practical search for these D° bound 
states. 



3 Results 

We look for D°-nucleus bound states by solving the Schrodinger equation: 



V 2 

-- + K pt (r) 



# = (B-iT/2)^, 



(6) 



where B is the binding energy (B < 0), T/2 the half- width of the bound state and p is the D- 
nucleus reduced mass. As mentioned above, we will present results from two different potentials 



opt l 



V„(r,E = B) 



and V opt (r) 



V D (r,E = 0). 



We solve the Schrodinger equation in coordinate space by using a numerical algorithm 
which has been extensively tested in similar problems of pionic [28| [3Tj and antikaonic [32] atomic 
states and in the search of possible antikaon [32] and r\ [23] nucleus bound states. Charge densities 
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Figure 2: Binding energies and widths for different -D°-nucleus states obtained using the strong energy 
dependent D— potential. 



are taken from Ref. [33]. For each nucleus, we take the neutron matter density approximately 
equal to the charge one, though we consider small changes, inspired by Hartree-Fock calculations 
with the DME (density-matrix expansion) [34J and corroborated by pionic atom analysis [35] . In 
Table 1 of Ref. [32] all the densities used throughout this work can be found. However, charge 
(neutron matter) densities do not correspond to proton (neutron) ones because of the finite size 
of the proton (neutron). We take that into account following the lines of Ref. [28J and use the 
proton (neutron) densities in our approach. 

Results, binding energies and widths, from the energy dependent potential are shown in Table [TJ 
These results are also presented in Fig. [5J where the states found for different nuclei and orbital 
angular momentum are collected. The bound states for which the absolute value of the binding 
energy is much smaller than the corresponding half width have not been presented, because they 
mix with the continuum energy spectrum and do not define clear bound states. For instance, in 
40 Ca we also find a 2s state with a binding energy of —1.3 MeV and half-width of 6.7 MeV, which 
mixes with the continuum. We have disregarded the consideration of this state and others alike. 

From the results of Table HJ we conclude there exist chances to see distinct peaks corresponding 
to D° bound states. Let us consider angular momentum I = 0. For medium size nuclei, up to 
40 Ca, there exists only one relevant D° bound state in each nucleus, with a binding energy around 
— 10 MeV and half- width roughly below 5 MeV, which can be subject to experimental detection. 
For heavier nuclei, like 118 Sn and 208 Pb, we find two I = bound states, with an energy separation 
similar to their half-widths (around 4 MeV). In the case of angular momentum 1 = 1, and for 
medium size nuclei (from 24 Mg up to 32 S) we find only the lp state, with \B\ smaller or equal 
than the T/2. Thus, these states are in the edge of the possible experimental determination. For 
heavier nuclei like 40 Ca and 118 Sn there is a well defined lp bound state. In the case of 208 Pb, we 
find the lp and 2p bound states separated by 4.3 MeV, while their half-widths are 3.6 and 4.9 
MeV, respectively . We also find states with I = 2 for 118 Sn and 208 Pb, and I = 3 for 208 Pb. 
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In summary, for light and medium nuclei (from A = 12 up to A = 40 ), there are observable 
D° mesic nuclei Is states. For heavier nuclei, and in addition to the Is level, we also find clearly 
observable lp states assuming that different angular momentum can be separated. For all cases, 
the bound D° meson is not orbiting around the nucleus, but rather it is embedded inside of it. For 
instance this can be seen in the left panel of Fig. EJ where the squared absolute values of the radial 
wave function for the Is and lp levels in 208 Pb, together with the nuclear density, are shown. 

To clarify the most relevant aspects of the model, we have also obtained the bound state 
spectrum when some of the ingredients of the full model are not considered. In Table |2j we show 
results for 12 C, 40 Ca and 208 Pb nuclei with different approximations. First, we examine the effects 
produced by the strong energy dependence of the optical potential and have computed energies 
and widths with the energy independent optical potential of Eq. (j2J). We see that when the energy 
dependence of the potential is neglected, the states become more bound by a factor ranging from 
1.6 in 12 C to 2 in 208 Pb, and also the widths of the states became larger by approximately a factor 
of two. These results can be easily understood by looking at the energy dependence of the optical 
potential in Fig. [TJ There, we see that both, the real (attractive) and imaginary (negative) parts 
of the optical potential decrease, in absolute value, by a factor larger than 2 from the threshold, 
q° = mo, to q° = mr> — 15 MeV. Thus, the energy dependence of the optical potential plays a 
major role and need to be taken into account. 

Next, we examine in Table |2] the importance of the selfconsistency in the calculation. To 
this end, we have recalculated the full spectrum by using an energy dependent optical potential 
deduced from the D selfenergy in the low density limit, Ili ow ~ Tp, where T is the ND T-matrix 
in free space. This amounts to use T p=0 , instead of T p in Eq ([1]). In Fig. [TJ we also compare the D 
selfenergy calculated for normal nuclear matter density with and without selfconsistency for zero 
momentum and energies below threshold. Though, the real part of the Yl£,(q° ,0, p) selfconsistent 
selfenergy does not differ much from that of the low density theorem approach, however the 
absolute values of the imaginary part of the II low approach are around a 60% larger than the exact 
result of Il£)(g ,<f = 0,p) for p = p , and in the energy region —20 MeV < q° — tud < 0. As 
a consequence, the obtained bound states (see Table [2]) using the low density limit have larger 
widths and smaller binding energies (due to the repulsive effect of the larger imaginary part) than 
those calculated with the exact selfconsistent potential. 

After analyzing the relevance of the ingredients of the model, we also compare our results in 
208 Pb with those obtained in Ref. [§]. There, a relativistic mean field calculation is carried out 
that leads to binding energies much larger than those obtained here. For instance, the Is state in 
208 Pb is bound by almost one hundred MeV and has no width in the model of Ref. [9], whereas we 
find a binding of only about ten MeV with a width of 6.6 MeV for the same Is level (see Table [2]). 
Of course, our coupled channel unitary and selfconsistent model is much more elaborated and it 
is also able to predict decay widths. 

Finally, we address the possibility of finding D + mesic nuclei. The calculated D— nucleus 
strong optical potential is the same for D° and D + mesons because of the isospin symmetry. The 
found full optical potential is not very deep. The strength of the attractive potential is, at most, 
20 MeV (see Fig. [JJ). A positive charged D + meson also feels the nuclear attraction of the nuclei. 
However, the Coulomb repulsion due to the positive electric charge of the nuclei is important, 
especially for heavy nuclei. We have considered a light nuclei, 12 C, and studied the Z) + -nucleus 
bound states by adding the Coulomb repulsive potential to the optical one. We have found that 
the Is state has a binding energy of B = —4.6 MeV and a half-width of T/2 = 6.4 MeV which is 
larger that \B\. Hence this state is not a clear case for experiments because it has a significant 
overlap with the energies of the continuous energy spectrum. Obviously for heavier nuclei, with 
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Table 1: {—B, T/2) in MeV for D°— nucleus bound states calculated with the energy dependent 
selfenergy. 





12 C 


24 Mg 


27 A1 


28 Si 


32g 


40 Ca 


118 Sn 


208p b 


Is 


(7.0, 5.0) 


(8.4, 4.1) 


(10.3, 5.1) 


(10.3, 5.0) 


(10.3, 4.7) 


(10.7, 4.5) 


(11.3, 3.6) 


(11.5, 3.3) 


Ip 




(3.5, 4.9) 


( 5.0, 6.5) 


( 5.3, 6.4) 


( 5.8, 5.9) 


( 6.8, 5.4) 


( 9.5, 4.0) 


(10.2, 3.6) 


Id 














( 7.2, 4.7) 


( 8.7, 4.0) 


2s 














( 6.4, 4.8) 


( 8.1,4.2) 


If 
















( 6.9, 4.6) 


2p 
















( 5.9, 4.9) 







Table 2: (-B, T/2) in MeV for D° -nucleus bound states in 12 C, 40 Ca and 208 Pb. The first 
set of results have been obtained by using the energy independent potential of Eq. (P2J). In the 
second set of results, though the energy dependence of the optical potential is kept, binding 
energies and widths have been obtained from a low density optical potential, where effects due 
to the selfconsistency treatment have been neglected. The last column gives, for comparison, the 
predicted binding energies of Ref. j9] for 208 Pb (in MeV). 





Energy independent V op t 


Energy dependent V op t 


without selfconsistency 


-B of Ref. [9] 




12 C 40 Ca 


208p b 


12 C 40 Ca 


208p b 


208p b 


Is 


(11.1, 10.9) (18.4, 13.5) 


(20.2, 11.0) 


(7.5, 10.2) (12.2, 7.5) 


(13.0, 5.2) 


96.2 


Ip 


(11.0, 11.4) 


(17.9, 10.7) 


(7.3, 11.1) 


(6.6, 9.5) 


93.0 


Id 




(15.0,10.3) 




(2.6, 12.8) 




2s 




(13.9,10.1) 




(1.5, 13.7) 


88.5 


If 




(11.7,9.8 ) 








2p 




(9.9,9.5) 









more protons, there is not a better chance of finding D + — nucleus bound states. We do not study 
the case of nuclei lighter than 12 C, because the LDA used here is not reliable enough for small 
nuclei, specially when significant cancellations among Coulomb and optical potential are taking 
place. Summarizing, with the optical potential considered in this work, there are no chances of 
finding D + mesic nuclei with atomic number Z equal or larger than 6. 

In the coupled channel calculation of Ref. [25], used as input for the D meson optical potential 
of this work, the nuclear matter vector meson D* selfenergy was calculated as well. The in-medium 
D* selfenergy is found there to be repulsive for energies around the threshold and densities p < po, 
hence that model predicts that there are not _D*-nucleus bound states. 

The SU(8) model used here as a kernel in the coupled channel unitary calculation reduces to the 
SU(4) model for J = 1/2, when the vector meson coupling constants are set to zero. To establish 
the relevance of the inclusion of the vector mesons, specially of the D*, in the D meson dynamics, 
we have also performed the calculations with the SU(4) model. In the left panel of Fig. HJ both 
the SU(8) and SU(4) D selfenergies are displayed for comparison. There, it can be appreciated 
that at threshold (q° = mr>) and for p < p , the SU(4) potential is small and repulsive, while the 
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Figure 3: Left: Squared absolute value of the radial part of the wave function of the Is and lp levels 
in 208 Pb. The nuclear density for lead is also shown. Right: Possible mechanism for production of D° 
mesic nuclei using a beam of antiprotons. 



SU(8) model provides a small attraction. However, for q° — vtld = —15 MeV and p < p , both 
potentials are attractive. However, the imaginary parts of both selfenergies are quite different, 
being much larger for the SU(4) case (about —60 MeV for p = po) than for the SU(8) model. This 
is due to the behaviour at finite density of the different resonance-hole structures of the D-meson 
selfenergy close to the DN threshold in the SU(4) [TTJ [19] and SU(8) [25] models. The large 
imaginary part of the SU(4) optical potential induces an effective repulsion, which together with 
the small attraction of its real part leads to the no existence of D°— nucleus bound states within 
the SU(4) model. This is in contrast to the SU(8) model that turns out to be attractive enough 
to give the bound spectra shown in Table [TJ and Fig. [2J 

A word of caution must be said here. The above discussion overlooks the fact there will be 
some extra contributions to the D— selfenergy from the DN — > D*N interaction mediated by pion 
exchange from the p— wave D* — > Dn and NNtc vertices (see right panel of Fig. H]). At threshold, 
these extra terms will be purely real, since there is no phase space for the reaction DN — > D*N, 
and they would also affect to the s— wave free space amplitudes derived in Ref. [26]. Because of the 
p— wave nature of the involved vertices, we expect the s— wave part of these contributions to be 
small near threshold, and their contribution to be effectively accounted for by the renormalization 
parameters used in [26], that were adjusted to reproduced the lowest-lying s— wave charmed 
baryon resonances. Note, however that these new mechanisms will be source of extra imaginary 
part in the case of the D*— selfenergy. 



4 Conclusions 

We have used a recent D°— nucleus optical potential, evaluated within a self-consistent unitarized 
coupled-channel approach to find D°— nucleus bound states. We have shown that selfconsistency 
effects and the energy dependence of the optical potential play major roles, and need certainly 
to be taken into account. The potential is attractive and we find bound states for all studied 
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Figure 4: Left: Comparison of the selfenergy over 2m p versus nuclear density for different 
values of the meson energy q° = nip and q° = m^— 15 MeV for SU(8) and SU(4) models. Right: 
Contribution to the D— selfenergy induced by the p— wave D* — > Dn and NNir vertices. 



nuclei through the whole periodic table. On the other hand, it produces states with relatively 
small half widths, smaller in general than the binding energies, and in many cases smaller than 
the separation among levels, which makes possible their observation. 

The strong energy dependence of the potential is due to the large effect on the D selfenergy 
of the charmed resonances, and their medium modifications [25], appearing close to DN thresh- 
old [26]. Taking into account the energy dependence reduces the strength of both the real and 
imaginary parts of the potential, below the D° threshold. This leads to substantially narrower 
states, but at the same time also with smaller binding energies. The heavy nuclei accommodate 
many D°-bound states and the separation of the levels, for a fixed angular momentum, is about 
half width of the states. The best chances for observation of bound states are in the region of 
24 Mg, provided an orbital angular momentum separation can be done, where there is only one s— 
bound state and its half width is about a factor of two smaller than the binding energy. In any 
case we would like to stress out that even if no broad bumps are found in the experiments, they 
should find some strength in the bound region stretching in energy down to the sum of the binding 
energy plus the half width of the bound states. Short of having the values for the binding energy 
and width of the states, this more limited information is still very valuable to gain some knowledge 
on the D nucleus optical potential, and it should stimulate experiments in this direction. This can 
be of particular interest for the future PANDA@FAIR physics program. Nevertheless, we should 
point out that the production and the experimental detection of these D°— nucleus bound states 
will likely be a quite difficult task. One might think in reactions of the type (see right panel of 
Fig. E} 

p + A z D- + {(A-l) Z -D\ Qmid 
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to be investigated in the future facility FAIR. Since the above reactions are two-body— j-two-body, 
the outgoing D energy is fixed, for a certain scattering angle in the laboratory system, and the 
creation of the D° bound state will be identified as a peak in the da/dQdE(D) over a background 
of inclusive (p, D) cross section. The High- Energy Storage Ring (HESR) at FAIR running with full 
luminosity (limited by the production rate of 2 x 10 7 p/s) at momenta larger than 6.4 GeV would 
produce about 100/s D meson pairs around t/>(4040) [35]. Assuming that one per ten million^] of 
the produced D° mesons is trapped in a bound orbit, we are then left with a production rhythm 
of around only one event every 10 5 seconds or equivalently to few hundred events per year. 

Recoilless production reactions have proved to be more efficient in the case of detecting deeply 
pionic bound states [3E1 ESS HO]. One possible reaction, where negligible momentum transfers 
could be achieved, is 



From the theoretical point of view we would expect sizeable formation peaks over a flat background. 
However, since the D* is an unstable particle, it will not be possible to create a D* beam, which 
makes, in practice, unfeasible the above reaction^. The above discussion brings us to reconsider 
the reactions of Eq. (j7j), and creating in the primary pN collision, not only a DD pair, but also 
a virtual pion, which will induce many body modes (particle-hole) in the nuclear medium [12]. 
Such modes will carry almost no energy, but high momentum components, which would allow the 
virtual D° meson in the right panel of Fig. [3] to have a significantly smaller momentum, being in 
this way, significantly increased the D° bound state production cross section. The main drawback 
would be that one is not facing now with a two-body going to two-body reaction. Hence to 
guaranty the creation of the D°— nucleus bound state, it would be needed to look at the decay 
products, A c 7r and S c 7r pairs, after the absorption of the D° by the nucleus. However, once the 
decay products of D°— nucleus bound state need to be detected, one realizes that it is unnecessary 
to use antiprotons (secondary beam) as projectiles, and instead one can benefit from the use of 
protons (primary beam), obtaining in this way a large enhancement factor (~ 10 5 ) in the incoming 
flux. Thus, one could look at reactions of the type 



Given the fact that there are more that two particles in the final state, it is possible to pick 
up regions of the phase space where the momentum, that needs to be accommodated in the D° 

1 This fraction should be understood only as an educated guess. In reactions of the type of that of Fig. |31 
the momentum transferred q to the bound D° is fixed. The cross section for the reaction is then proportional 
to \& n i m (q)\ 2 , where ^ n im{q) is the D° bound wave function in momentum space |37[ 138] . Since in the reactions 
of Eq. ([7]) the momentum transfer is quite large, around 2 GeV at threshold, and at least 1 GeV for incoming 
antiproton momenta of 10 GeV or larger, one can immediately see that the cross sections will be small, since the 
D° bound wave function has difficulty to accommodate such large momenta; but the reaction is possible. Given 
the typical size of the D° bound states (see for instance the left panel of Fig. [3]), momentum transfers of about 
0.2 GeV can be easily accommodated, and in these circumstances one might expect that one per cent of the D° 
mesons could be trapped |37| 138], The extra six orders of magnitude of reduction, implicit in the fraction 10~ 7 , 
would account for the expected large reduction of the D° bound wave function when the momentum increases from 
0.2 GeV to 1 GeV (we have assumed a factor 10 -3 ). 

2 Around ten millions of 5Wfl'*' pairs per year are expected to be produced at BESIII in the region of the 
^(3770) and V(4160) resonances |41j . Even if it were possible to put a thin nuclear target of 1 cm 2 at a distance 
at small as 1 meter from the collision area, we would not have more than one per year charmed meson interacting 
with the nuclear target. 






(9) 



M> A C 7T + X 
S C 7T + X 
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bound wave function, would be sufficiently small to make significantly larger the probability of 
D° trapping. Moreover, the above reaction is coherent; there is no change of charge in p + iV — > 
p + D° + D° + iV and the final nuclear state is the same as the initial one. This provides a factor 
A 2 in the cross section versus a factor A in the inclusive reaction which gives the background [3H] • 
Because of that, D°— bound states in heavier nuclei than magnesium, mentioned above, might 
have better chances to be detected. 

Finally, we would like to point out that in this work we have used a model for the in medium 
D— meson selfenergy that, it is based on a SU(8) spin-flavor extension [26] of the SU(2) Weinberg- 
Tomozawa (WT) pion-nucleon s— wave interaction in the free spaced. When vector mesons are not 
considered, our in medium selfenergy reduces, up to minor details, to that deduced in Ref. [T9] . 
This latter one is based on a SU(4) flavor extension of the SU(2) vacuum WT pion-nucleon 
s— wave interaction. The optical potential deduced from Ref. [19] is not attractive enough to 
bind D° nuclear states. This is in sharp contrast with our findings here. One of the major 
differences between both approaches, and mostly responsible for this distinctive difference, is that 
the model used here treats heavy pseudoscalar and heavy vector mesons on equal footing, as 
required by Heavy Quark Symmetry (HQS). This latter symmetry is a proper QCD spin-flavor 
symmetry [431 HH H5] when the quark masses become much larger than the typical confinement 
scale, Aqcd- 
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